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In this work, we examine the finite temperature properties of the CPT-even and Lorentz- 
invariance-violating (LIV) electrodynamics of the standard model extension, represented by the 
term W a v PV F al ' F pip . We begin analyzing the hamiltonian structure following the Dirac's procedure 
for constrained systems and construct a well-defined and gauge invariant partition function in the 
functional integral formalism. Next, we specialize for the non-birefringent coefficients of the tensor 
Wa Vpv - In the sequel, the partition function is explicitly carried out for the parity-even sector of 
the tensor Wavptp- The modified partition function is a power of the Maxwell's partition function. 
It is observed that the LIV coefficients induce an anisotropy in the black body angular energy den- 
sity distribution. The Planck's radiation law, however, retains its frequency dependence and the 
O |. Stefan-Boltzmann law keeps the usual form, except for a change in the Stefan-Boltzmann constant 

by a factor containing the LIV contributions. 
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^ ; I- INTRODUCTION 

Q-: 

The researches about Lorentz and CPT violation are commonly performed under the framework of the standard 
model extension (SME) developed by Colladay and Kostelecky [lj. The SME is an enlarged version of the usual 
standard model that embraces all Lorentz-invariance-violating coefficients whose tensor contractions yield Lorentz 
scalars in the observer frame, and in the particle frame are seen as sets of independent numbers. A strong motivation 
to study the SME is the necessity to get some information about underlying physics at Planck scale where both 
the Lorentz and CPT symmetries can be broken due to quantum gravity effects such as it is suggested by string 
theory El. The photon sector of the SME has been intensively studied in the latest years with a double purpose: 
the determination of new electromagnetic effects induced by the LIV interactions and the imposition of stringent 
upper bounds for the magnitudes of the LIV coefficients. Such investigations have connections with the Carroll- 
Field- Jackiw electrodynamics Ufij], consistency aspects 0], polarization deviations for light traveling over large 
cosmological distances El, El @] , Cerenkov radiation [8| , radiative corrections Q , electromagnetostatics and classical 
solutions 0] , U | , 0] , [l 3l| . radiation spectrum of the electromagnetic field and CMB 
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and quantum el^dynlmfos processes Q, 0,0, 0,0,0, and synchrotron radiation Q- For a large and 



interesting review on the photon sector and related issues, see Ref. 23]. Lorentz violation and its implications have 
been studied in several diverse respects 0],0] and also in other theoretical environments 0|. 

The most general renormalizable form of the Lorentz-covariance- violating electrodynamics of the SME photon sector 
can be expressed by the following Lagrangian density 

C = --F av F av - (k AF ) A a F pip - ~W av ™ 'F av F pv , (1) 

where e^ apip is the totally antisymmetric Levi-Civita tensor (e 0123 = 1), F av is the electromagnetic field tensor, A a is 
the vector potential, [kAF)g — (0, k^) has the dimension of mass and describes a super-renormalizable (dimension 
3) coupling, W avpv is a renormalizable, dimensionless coupling giving raise to a dimension 4 operator. The tensor 
W avpip has the same symmetries of the Riemann tensor \W avpip = —W vapip ,W aupip — —W avipp , W avpip = W pipav \ and 
a double null trace which yields only 19 independent components. 

The term e^ KX (k AF )^ A V F^\ is CPT-odd and it was first introduced by Carroll-Field-Jackiw (CJF) [3], who studied 
the modifications produced by this term (including vacuum birefringence) in the classical Maxwell electrodynamics. 
It yields a causal, stable, and unitary electrodynamics only for a purely spacelike background (see work of Adam & 
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Klinkhamer in R ef. [ Hi)- The term W fJ,VK,x F pu F K \ 1 composed of 19 elements, is CPT-even and was much investigated 
in Refs. 12, 13, 17, 20, 21, 2^, 27|. It yields an electrodynamics not plagued with stability illness. 

As the LIV terms alter the light propagation, it is natural to infer that the thermodynamical properties of the theory 
is modified as well. In a recent work 14j . it was investigated the influence of the CFJ term on the thermodynamics 
of the Maxwell field, using the usual formalism of finite temperature field theory [jjj . It was first analyzed the 
Hamiltonian structure of the model using the Dirac formalism in order to define the partition function of this theory 
without ambiguities. In the sequel, the LIV corrections induced on the black body Planck distribution were carried 
out, properly examined, and related to the cosmic background radiation (CMB). A similar investigation remains to 
be done for the CPT-even photonic sector of the SME, being it the main purpose of the present work. 

This paper is organized as follows. In Sec. II, we discuss the Hamiltonian structure of the CPT-even electrody- 
namics, investigating the constraints structure of the theory (by means of the Dirac formalism). In Sec. Ill, we write 
the general partition function of this electrodynamics into the functional integral formalism by using the constraint 
structure and the gauge fixing conditions established in Sec. II. This partition function is then particularized and 
explicitly evaluated for the parity-even sector of the theory. The modified partition function is written as a function 
of the usual Maxwell partition function. In Sec. IV, we present our final remarks discussing the modifications induced 
on the Maxwell theory and comparing it with the Carroll-Field-Jackiw model at finite temperature results. In the 
Appendix, we evaluate the dispersion relations for the parity-even sector of the CPT-even electrodynamics, which 
corroborate the results obtained for the partition function of this work. 

II. THE CPT-EVEN AND LIV ELECTRODYNAMICS OF THE STANDARD MODEL EXTENSION 

In the present work, we just study the CPT-even and LIV electrodynamics of the SME, so that we will consider 
(kAF)p = 0. Therefore, the Lagrangian density given by Eq. is reduced to 

H = -\p av F av -\w avpv F av F pv , (2) 
which yields the following Euler-Lagrange equation for the gauge field 

d v F vp - W» vp<p d v F ptp = 0. (3) 

A. The Hamiltonian structure 

In order to accomplish the Hamiltonian analysis of this model, we begin defining the canonical conjugate momentum 
of the gauge field as 

with which we can write the fundamental Poisson brackets (PB): {A p (x) , n v (y)} — S p 5 (x — y). 

From the Eq.Q, it is easy to note that tt° — 0. Such a null momentum yields a primary constraint <fii — ir° w 
(into the Dirac formalism, the symbol ~ denotes a weak equality). Also, the momenta n k arc defined via the following 
dynamic relation 

n k = D kj F oj -W ok ^ l F jl , (5) 
where the nonsingular and symmetric matrix D k j is defined by 

D k j = Skj - 2W feoj. (6) 

Then, the velocities A k are given as 

A k = d k A + (D- l ) k . [tt* + W^ mn F mn ] , (7) 
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while the canonical Hamiltonian density is explicitly written as 

H c = \ [n k + W okmn F mn ] (D~ l ) [y + W 0jmn F mn ] + n k d k A + \ {F jk f + \w h ? lm F kj F lm . (8) 



Following the usual Dirac procedure, we introduce the primary Hamiltonian (Hp) by adding to the canonical Hamil- 
tonian all the primary constraints, Hp = He + J d 3 y Cn° , where C is a bosonic Lagrange multiplier. The consistency 
condition of the primary constraint, tt° — {ft , Hp} « 0, gives a secondary constraint 



<2 



= d k n k «0. (9) 



It means that the Gauss's law structure is not modified by the CPT-even and LIV background. Nevertheless, express- 
ing it in terms of the electric and magnetic fields, we can note the explicit coupling between the electric and magnetic 



sectors even in the electrostatic regime 

The consistency condition of the Gauss's law gives <j>2 = 0. Thus, the secondary constraint is automatically 
conserved and there are no more constraints in this model. The bosonic multiplier of the primary constraint remains 
undetermined, being an evidence for the existence of first-class constraints. This is verified by computing the PB 
between the primary and the secondary constraints: {7r°, dkft k } = 0. The constraints (j>\ = tt° ~ and 02 = d k ft k ~ 
reveal that the CPT-even and LIV electrodynamics has a similar constraint structure as the Maxwell electrodynamics. 

B. Equations of motion and gauge fixing conditions 

Following the Dirac conjecture, we define the extended Hamiltonian (He) by adding all the first-class constraint to 
the primary Hamiltonian, 

H E = H c + J dy [Cfa + A0 2 ] ■ (10) 

Under this Hamiltonian, we compute the time evolution of the field variables of the system 

A a = {A ,H E } = C, (11) 

A k = {A k ,H E } = (D- 1 )^ [it* + W°' mn F mn ] + d k A - d k A, (12) 

Both equations show that the dynamic of gauge field A^ remains arbitrary. However, the second equation is similar 
to the Lagrangian equation ([7]) if and only if A = 0. Thus, we should impose a gauge condition in such a way to 
fix A — 0. As it is well-known, the Dirac algorithm requires a number of gauge conditions equal to the number of 
first-class constraints in the theory. However, those gauge conditions must be compatible with the Euler-Lagrange 
equations, such that they should fix A = and determine the Lagrangian multiplier C. The gauge conditions together 
with the first-class constraints should form a second-class set. 
From the equation of motion for Aq, 

D jk djd k A - II . ,,,;.<), /•;,;. - d (DjkdjAk) = 0, (13) 

we set as our two gauge fixing conditions 

V>i = DjkdjAk m 0, i> 2 = DjkdjdkAo - WoijkdiFjk r> 0. (14) 

The consistency condition for ifji gives Dj k djd k A — 0, which fixes A = 0. The consistency condition for ip 2 gives 
an equation for the multiplier C 

D jk djd k C - Wo^kdtF^ » . (15) 

Consequently, we have determined all the Lagrange multipliers. Therefore, the set S a = {4>i, 4>2, ipi, ^2} is a 
second-class one. 
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The next step is to compute the Dirac brackets to know the field algebra. Thus, after a long computation, we find 
that the non-null Dirac brackets are 

{A k {x) , n j (y)} D = S k] 6 (x - y) + F> m jcf m d%G (x - y) , (16) 

{A (x),n k (y)} D = -2W 0ijk d?d?G{x-y), (17) 

where the matrix Dj k is given by Eq. ([6]) and G (x — y) is the Green function for the modified Poisson equation: 

D jk d j d k G(x.-y) = -8(x-y). (18) 

The Dirac brackets for the physical degree of freedom (fl6|) do not reflect the transverse character of the gauge field, 
however if we choose d k A k ~ as a gauge condition, the DB is reduced to the usual transverse commutation relation 

(Qx \ 
J^(x-y). (19) 

Here, we need to observe that, at quantum level, the transverse character of the gauge field can be explicitly proven by 
computing the Ward identity for its 1PI 2-point function, W {x - y) = (Ug^ - d^d" + ^d^d" - S^) S {x - y), 
thus d^ v [x-y) = ^ 1 d v 5 (x-y). 

Under the Dirac brackets, the canonical Hamiltonian ([5]) reads as 

H = J dy S^E k D k jE j + ^B 2 + ±W kjlm F kj F lm } . (20) 

In general, for a sufficiently small W^ 1 " 7 , a positive-definite Hamiltonian is guaranteed, thus providing a stable 
quantum theory and a well-defined partition function associated with the CPT-even and LIV electrodynamics. Now, 
we proceed to the computation of the partition function, performing the analysis of its implications to the black body 
radiation problem. 



III. THE PARTITION FUNCTION 

The next step is to study the thermodynamical properties of the CPT-even photon sector of the SME. The funda- 
mental object for this analysis is the partition function. The Hamiltonian analysis performed in the previous section 
allows to define the partition function (in a correct way) into the functional integral representation 

Z{(3) = JvA^S^S^S^S^) |det{E a (a;),S 6 (y)}| 1/2 exp|y& {m*d T A„ -H C )J, (21) 

where E a = {4>\, (j)2, ipi, "02 } is a second-class set formed by the first-class constraints and the gauge fixing conditions, 
M ab {x,y) = {S a (x) , T, b (y)} is the constraint matrix whose determinant is det (— Djkdjdk) ■ Given the bosonic 
character of the gauge field, its functional integration can be performed over all the fields satisfying periodic boundary 

conditions in the t— variable: A^ (r, x) = A^ (r + /3, x). The short notation / dx denotes / dr d 3 x, and Tie is the 

J p Jo J 

canonical Hamiltonian given by Eq. ([8]), and [3 — 1/ksT, where Ub is the Boltzmann constant. 

By performing the integrations over the canonical conjugate momenta and doing the following redefinitions: F rk — 
d T A k - d k A T = -F kT and 

WokOj — —W TkT j , W() krnn — iW T kmn, (22) 

we find the partition function for the CPT-even photonic sector of the SME as 

Z {(3) = TV det {-Djkdjdk) J VA a 5 {D jk djA k ) exp I^Jdx - ^F ab F ab - ~ W abcd F ab F cd ^ , (23) 
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where a,b,c,d = t, 1,2,3. This partition function is not explicitly covariant. However, it is well-known that if the 
covariance is explicit, the calculation process becomes more manageable. The procedure to pass from a non-covariant 
gauge to a covariant one (like the Lorentz gauge d a A a = 0) can be performed using the Faddeev-Popov ansatz. Thus, 
choosing the following Lorentz gauge G [A a ] = — £,~ x ^d a A a + /, where / is an arbitrary scalar function and £ is a 
gauge parameter. Therefore, after some algebra, we find the partition function to be 



Ztf) = J DA a det (-j^j exp dx - l -A a -Ub ah - Q - 1^ d a d b + S ab A b | 



(24) 



where □ = d a d a — (d T ) + V 2 . We have also defined the symmetric LIV operator 

Sab = 2Wacdbd c d d = S ba . (25) 

For convenience, we choose the Feynman gauge £ = 1. Performing the gauge field integration, we find 

Z 09) = det (-□) [det {-nS ab + S ab )]~ 1/2 • (26) 

It is illustrative to mention that this result is similar to that obtained for the Carroll-Field-Jackiw electrodynamics, 
where S ab = e acdb (K AF ) c d d . 

Given the high complexity of the CPT-even term, in order to turn feasible the explicit evaluation of the partition 
function, the tensor W acdb should be specialized for simpler configurations. It is done at zero temperature in Refs. 
@ 5 0|j fr° m which one knows some useful parametrization for the tensor Wuvafi i n terms of four 3x3 matrices, 
kde, k hb, k db, khe : 

( KDE y k = -2W°i° k , (K HB y k = l e 3P q E klm W pqlm , (K DB ) jk = - (K HE ) kj = e kpq W°i pq . (27) 

The matrices kde and khb contain together 11 independent components, while kdb and khe possess together 8 
components, which sums the 19 independent elements of the tensor W acdb - Such coefficients can be parameterized in 
terms of four traceless matrices and one trace element. The parity-odd sector is written as 

(«o+)fcj = ~(kD-B + KHE)kj, («o-)fcj = ^DB - Kff E )kj, (28) 

while the parity-even sector is read in terms of two matrices and one trace element, 

(Ke+)kj = 7;( K DE + K HB )kj, («e-)fcj = - K HB )kj - nS k j , n = ^tr (kde) ■ (29) 

The matrix k + is antisymmetric while the other three are symmetric. Ten of the 19 elements of the tensor W avp ip 
(5 belonging to k q _ and 5 to K e +) are strongly constrained by birefringence data (at the level of 1 part in 10 32 ) 
@, 0, [23] • From the nine remaining nonbirefringent coefficients, three are contained in the parity-odd matrix k +- 
The parity-even sector encloses six elements (five in the matrix K e _ and the trace element, n). 
The prescriptions (|2T|) . taking into account the finite temperature redefinitions (|2"2"|) , are read as 

i K DE) k j — 2W T k T j, {nHB) k j = ^kpq e jmnW pqmn , (KDB)ki = ~ ( K HE)j k = W T kpq£jpq- (30) 

We should now carry out the determinant of the operator (—\D5 ab + S ab ) for the six non-birefringent components 
of the parity-even part of the W acdb tensor. 



A. The parity-even sector 



The parity-even sector is composed of an isotropic component and five anisotropic components - the elements of 
matrix /? e _. We now evaluate the partition function for this sector. 
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1. The isotropic contribution 



We first isolate the isotropic part of the parity-even sector by imposing {n e -)j k = 0, retaining only the component 
n. The functional determinant for the operator (— D6 a b + S a b) is now given as 



det (-nSab + S ab ) = det (n + if [-D] 2 det 



2n , 
n+ 1 



while the partition function becomes 



-Trln 



2n 

i 

71 + 1 



We can evaluate the involved trace by writing the gauge field in terms of a Fourier expansion, 





A a (T,*) = (~) $>^+ x 'P)i (n,p), 



(31) 



(32) 



(33) 



n,p 



2mr 



where V designates the system volume and uj n are the bosonic Matsubara's frequencies, ui n — , for n = 0, 1, 2, 
The contributions of the two modes of the gauge field are expressed as 

,2 l-n 



\ / fryi f^s-. 



1 + n 



(34) 



Here, it should hold \n\ < 1 for yielding a well-defined partition function. By performing the rescalingp,; — > p t 



we obtain 



1+71 
1 — 71 



3/2 



IsxZa, 



where Za is the partition function of the Maxwell's electrodynamics, given by 

\nZ A = =■ / oIujlu 2 In (1 - e-^) = V— 



45/3 3 



1 + n 
~ n' 



(35) 



(36) 



From p5p . we see that the LIV partition function is obviously a power of Za, 

Z (f3) = (Z A ) a(n) , (37) 

for a (n) = ((1 + n)/(l - n)f /2 . With this result, it is easy to show that both the modified Planck's radiation and 
the Stefan-Boltzmann's law of the isotropic sector are those of the Maxwell electrodynamics multiplied by the factor 
a (n). Here, the energy density per solid-angle element remains isotropic. 



2. The anisotropic contribution 

The anisotropic coefficients of the parity-even sector are represented by the terms of the matrix (/t e -) • They can 
be isolated by setting n = 0. For evaluating the functional determinant, we should express the matrix (/t e -) " m a 
suitable way. As the matrix (f? e _) is symmetric and traceless, it can be parameterized in terms of two orthogonal 3D 
vectors, a and b, as 

(«e-) ife = T ( a J b k + b J a k) , (38) 
with a • b = and det (K e -) = 0. Then, the functional determinant of the operator (— OS a b + S a b) is 

det (-DS ab + Sat) = det M - W) det (-D) 2 det (-□ - V + ) det (-□ - V 2 .) , (39) 
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where the operators and V 2 . are given as 



2 _ 4 (a ■ V) (b ■ V) + b 2 (a ■ V)" + a 2 (b ■ V)- 



+ 

72 



a 2 b 2 



VI = (a-V)(b-V). 
With all these definitions, the partition function becomes 



\nZ{/3) = -ilndet [-□ - V+] - ~ bidet [-□ - V 2 J\ 



(40) 
(41) 

(42) 



representing the contributions of the two polarization modes of the gauge field. Let us observe that if we consider 
only the first order contribution of the LIV background, we have 



V 2 h «(a-V)(b-V), V 2 _ = (a- V) (b • V) 



(43) 



It means that the dispersion relation at first order are the same for both modes of the gauge field, once both modes 
give the same contribution to the partition function at first order. For an alternative evaluation of the dispersion 
relations, see Appendix. This result is compatible with the statements of Ref. 23 1. 

Again, the functional trace is carried out by means of the Fourier expansion (f33|) of the gauge field. The contributions 
of the two modes of the gauge field are expressed as 



lnZ + ((3) 



M 3p ^ i « 2 



lnZ_(/J) = --V 



(2tt) 3 
(2tt) 3 



m— — oo 



a a b 2 (a.pr +a 2 (b. P r+4(a.p)(b.p) 
(Um ) + P H 



[4 - a 2 b 2 ] 



ln/3 2 [( Wm ) 2 + p 2 + (a-p)(b.p) 



(44) 



(45) 



m— — oo 



In order to perform the momentum integrations, we consider the following coordinate system: the vector a is aligned 
with the x— axis, the vector b with the y— axis, so that a x b points along the z— axis. Expressing the momentum 
in spherical coordinates [p = u) (sin cos 0, sin sin 0, cos 0)], we achieve the following mode contributions: 



lnZ+(/J) 



lnZ_ (/?) 



1 V 



, , , du u? V ln/3 2 



(io m y + Lu z i + 



a 2 6 2 sin 2 9 + 2ab sin 2 9 sin 20 



(w m ) z + u}' 2 ( 1 + -afesin 2 #sin20 



(46) 



(47) 



where c?£l = sin 9 d9d<fr is the solid-angle element, a = |a| and b = |b|. By doing a rescaling of the variable w and 
performing the summation, we obtain 

-3/2 



InZ+09) 
lnZ_ (f3) 



1 

8^ 

1 

8^ 



\nZ A dn [1 



\ 2 sin 2 9 + 2X sin 2 9 sin 20 



4- A 2 



1 



In Z A (ffl 1+ - A sin sin 20 



-3/2 



(48) 
(49) 



where we have defined A = ab. Taking into account the outcome of Eq. (|36|) . and noting that the angular integrations 
can be exactly solved, the partition functions become 

1 



lnZ+(P) = ±(4 -A 2 ) 1/2 lnZ A ((3), 



lnZ_(/3) = (4 -A 2 ) 



-1/2 



lnZ A ((3). 



(50) 
(51) 



For having a well-defined partition function, the product ab = A must be bounded as < ab < 2. Remembering that 
InZ (P) = \nZ + ((3) + \nZ_ {(3) ,we can also show that the LIV partition function can be written as a power of the 
Maxwell's one, 

z (/?) = (z A ) sw , 
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with 



1 



5(A) = ±(4-A 2 ) 1/2 + (4-A 2 ) 



2N-1/2 



(52) 



Similarly to the parity-odd case, we observe that the modified Planck's radiation law and the Stefan-Boltzmann's 
law are those of the Maxwell electrodynamics multiplied by the factor S (A). However, the energy density distribution 
per solid angle, 

-3/2" 



120/5 4 



1 + ^Asin 2 6»sin20 



-3/2 



1 



y sin z 9 + 2A sin z 9 sin 26 
4~ A 2 " 



(53) 



possesses an explicitly dependence on 6 and 9 which reveals a higher degree of the anisotropy induced by LIV, as it 
can be shown at leading order 



120/3 4 



2 - - A sin 26 sin 2 9 
2 Y 



(54) 



The A linear dependence of the energy density may lead to an attainment of upper-bounds on the K e _ parameters 
using polarization data of the cosmic microwave background. 



IV. CONCLUSIONS AND REMARKS 



We have initiated this work establishing the Hamiltonian structure of the CPT-even sector of the electrodynamics 
of the SME. The constraint analysis allows the construction of a well-defined partition function which is given in 
(1261) for an arbitrary and sufficiently small tensor W a bcd- At once, we specialize our analysis for the non-birefringent 
components of the parity-even parts, for which we compute exactly the partition function. The expression (|52j) shows 
that it is a power of the partition function of the Maxwell electrodynamics, where the power is a pure function 
the LIV parameters. This way, the Planck radiation law retains its usual functional dependence in the frequency 
and the Stefan-Boltzmann law remains the same one, apart from a multiplicative global factor containing the LIV 
coefficients. It is observed that the LIV induces an anisotropic angular distribution for the black body energy density 
for the anisotropic parity-even (rt e -) coefficients. The anisotropic character of the angular radiation distribution 
reflects local energy density variations in relation to the Maxwell pattern induced by Lorentz violation. Despite such 
differences, the Stefan-Boltzmann law keeps the usual temperature behavior. This means that, notwithstanding small 
local fluctuations, the global radiation law maintains the T 4 — behavior. 

Since the LIV coefficients are constrained by very stringent upper bounds, the lower order non-null LIV contribution 
for the Maxwell thermodynamics would give a good information about the thermodynamical properties of the non- 
birefringent sector the model. It is observed that the isotropic contribution gives a linear correction in n, whereas 
the anisotropic contribution coming from the matrix H e _ only is manifest at fourth order, as it is shown by Eq. (1 5 [) . 
Hence, the pure anisotropic contribution is irrelevant when compared with the isotropic one. 

Moreover, we must highlight the differences between the thermodynamical properties of the CPT-even and the 
CPT-odd electrodynamics, first investigated in Ref. J3. Such difference stems from the Dirac's algebra of the 
physical variables. For the CPT-odd electrodynamics [14|, in the Coulomb gauge, we have attained 



{A k (x),TTj {y)} 



D 



4 



dkd-j 



{n k (x),^ (y)} D = - 



(k 



AF)i ■ 



^ {kAF) d lK 



v 2 



<5(x-y). 



(55) 
(56) 



Nevertheless, for the CPT-even case and for the Maxwell electrodynamics, the Dirac algebra is given only by Eq. 
(|55|) . The noncommutativity of the physical momenta, expressed in Eq. (I56p . is the fundamental reason for this sector 
to have different thermodynamical properties when it is compared with its CPT-even counterpart. Also, since the 
background Uaf is a dimensional parameter, the temperature dependence of the logarithm of the partition function at 
order (kAF) 2n changes as T 3 ~ 2n . It has as a consequence that the CPT-odd partition function can not be expressed 
as a power of the Maxwell one such as it happens in the CPT-even case. 
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A. APPENDIX: Dispersion relations 

In this Appendix, we write the dispersion relations for this CPT-even electrodynamics as a procedure to confirm 
the evaluation of the partition functions. A general evaluation for the dispersion relations may be developed from 
Eq. ([2|) and the matrix prescriptions (|27p . In terms of the matrices kde, k-db, k-hb the non-homogenous Maxwell 
equations (in the absence of sources) are 

djEj + (K DE ) ja d,E a - {n DB ) ja d,B a = 0, (57) 
d E k - ekjadjBa + {n DE ) kj d Ej - (n DB ) kb d B b - (n DB ) ab e bkj djE a - (n HB ) ab e bkj djB a = , (58) 

while the homogenous ones remain the same, doB k + e kab d a E b = 0, d a B a = 0. The wave equation for the electric field 
is 

(d t ) 2 E k - V 2 E k + tr (k DE ) ^ 2 E k + (K HB ) ab d a 8 b E k + d k 8 a E a + (n DE ) ka (8 t ) 2 E a 
-tr (kde) d k d a E a - (n HB ) kb d b d a E a + (n H B) ka V 2 E a - [nHB) ba d b d k E a 

+ (KDB) kc e c bad t d b E a - (n DB ) ab e bkc d c d t E a = 0. (59) 

We now specialize the wave equation (f59|) for the parity-even case, setting kdb = 0. We then express the matrices 
kde, k B b in terms of the K e+ and K e - 

(nDE) ab = (Ke+) ab + (Ke-) ab +n5a b , (60) 
[nHB) ab = (K e +) ab - (Ke-) ab -n5 ab . (61) 

Birefringence data impose (re e +) a6 = 0, so that 

(dt) 2 E k + n (dt) 2 E k + [Ke-) ka (d t f E a - V 2 E k + nV 2 E k - (K e ^) ka V 2 E a 

- (Ke~) ab d a d b E k + d k d a E a - nd k d a E a + (« e -) fc6 d b d a E a + (K e _) ba d b d k E a = (62) 
Retaining only the isotropic component (n 7^ 0,/5 e _ =0), we have 

(1 + n) (d t ) 2 - (1 - n) V 2 1 E k + (1 - n) d k 8 a E a = 0. (63) 



Using now the first Maxwell equation, (1 + n) d a E a — 0, we obtain (1 + n) (dt) — (1 — n) V z 
Fourier space, is reads as 



E k = 0, which in 



[(l + n)pg-(l-n)p 2 ]£ fc = 0. (64) 
This equation yields the following dispersion relation: 

(l + n)p3-(l-n)p 2 = 0. (65) 
This is the same expression contained in Eqs. (|31[ 132134}) . confirming our previous result: 



w± = ±|p|V(l-n)/(l + n). (66) 

Here, we see that the phase velocity associated with the modes of the photon field is the same, showing explicitly 
the nonbirefringent character of the isotropic coefficient of the parity-even sector, which is in full accordance with the 
statements of Ref. [23j. Moreover, we note the existence of positive and negative frequencies, uj + and u>-. 

We should finally consider the anisotropic components of the parity-even sector (n = 0,/5 e _ 7^ 0). In this case, the 
wave equation (|59|) reads as 



[□4a + (K e -) ka n]E a - [(K e -) cb d c d b 5 ka - (K e ^) kb d b d a ] E a = 0. (67) 
In momentum space, we have 

{ [p 2 - (Ke-) cb PcPb\ 5 k a + («e-)fcj, [p 2 Sab + PhPa] } E a = 0. (68) 
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We now use the same parametrization of Eq. (|38|) . where a and b are two orthogonal 3D vectors. Then, we have 



( K e-)cbPcPb = (a-p)(b-p), 

(K e _) fec [p 2 S jc +p c pj] = - {a k b } + ajb k )p 2 + -a k p 3 (b • p) + -b k pj (a • p) 



With it, Eq. d68j) is read as 



with 



M kj Ej = 0, 



Mkj = [p 2 - (a • p) (b • p)] 6 kj + - (a k bj + djb k ) p 2 + ^a k p 3 (b • p) + -b k pj (a • p) 
The dispersion relations are obtained from detM = 0. Computing the determinant, we get 



detM =( 1 - -a 2 b 2 ) p 2 [p 2 - (a • p) (b ■ p) 



4 (a ■ p) (b ■ p) + (a ■ p) 2 b 2 + a 2 (b ■ p) 2 
4 - a 2 b 2 



(69) 
(70) 

(71) 
(72) 

(73) 



By this way, we attain the exact dispersion relations 



P 2 = (a • p) (b • p) , 

p 2 = [4(a-p)(b-p) + (a-p) 2 b 2 +a 2 (b-p) 2 ][4-a 2 b 2 ]- 1 . 



(74) 
(75) 



These expressions confirm that the partition function for this case is the one stated in Eqs. (|42[ 144145)) . At leading 
order, these dispersion relations are the same one, 



w = ±lp| 



1 , l (a-p)(b-p) 
+ 2 Inl 2 



(76) 



implying absence of birefringence at leading order such as stabilished in Ref . [23j . 



Acknowledgments 



R. C. thanks Conselho Nacional de Desenvolvimento Cientffico e Tecnologico (CNPq) for partial support, M. M. 
F. is grateful CNPq and to FAPEMA (Fundacao de Amparo a Pesquisa do Estado do Maranhao) for partial support. 
J. S. R. thanks FAPEMA for full support. 



[1] D. Colladay and V. A. Kostelecky, Phys. Rev. D 55, 6760 (1997); D. Colladay and V. A. Kostelecky, Phys. Rev. D 58, 
116002 (1998). 

[2] V. A. Kostelecky and S. Samuel, Phys. Rev. Lett. 63, 224 (1989); 66, 1811 (1991); Phys. Rev. D 39, 683 (1989); 40, 1886 
(1989); V. A. Kostelecky and R. Potting, Nucl. Phys. B 359, 545 (1991); Phys. Lett. B 381, 89 (1996); V. A. Kostelecky 
and R. Potting, Phys. Rev. D 51, 3923 (1995). 

[3] S. M. Carroll, G. B. Field and R. Jackiw, Phys. Rev. D 41, 1231 (1990). 

[4] AA. Andrianov and R. Soldati, Phys. Rev. D 51, 5961 (1995); Phys. Lett. B 435, 449 (1998); A. A. Andrianov, R. Soldati 

and L. Sorbo, Phys. Rev. D 59, 025002 (1998); M. Frank and I. Turan, Phys. Rev. D 74, 033016 (2006); A. Hariton and 

R. Lehnert, Phys. Lett. A 367, 11 (2007). 
[5] C. Adam and F. R. Klinkhamer, Nucl. Phys. B 607, 247 (2001); Phys. Lett. B 513, 245 (2001); H. Belich, M. M. Ferreira 

Jr, J.A. Helayel-Neto, M. T. D. Orlando, Phys. Rev. D 67, 125011 (2003); -ibid, Phys. Rev. D 69, 109903 (E) (2004). 
[6] V. A. Kostelecky and M. Mewes, Phys. Rev. Lett. 87, 251304 (2001); V. A. Kostelecky and M. Mewes, Phys. Rev. D 66, 

056005 (2002). 

[7] V. A. Kostelecky and M. Mewes, Phys. Rev. Lett. 97, 140401 (2006); Astrophys. J. Lett. 689, LI (2008). 



11 



[8] R. Lehnert and R. Potting, Phys. Rev. Lett. 93, 110402 (2004); Phys. Rev. D 70, 125010 (2004); 129906(E) (2004); B. 
Altschul, Phys. Rev. D 75, 105003 (2007); Phys. Rev. Lett. 98, 041603 (2007); B. Altschul, Nucl. Phys. B 796, 262 (2008); 

C. Kaufhold and F.R. Klinkhamer, Nucl. Phys. B 734, 1 (2006); C. Kaufhold and F.R. Klinkhamer, Phys. Rev. D 76, 
025024 (2007). 

[9] R. Jackiw and V. A. Kostelecky, Phys. Rev. Lett. 82, 3572 (1999); J. M. Chung and B. K. Chung Phys. Rev. D 63, 105015 
(2001); J.M. Chung, Phys. Rev. D 60, 127901 (1999); G. Bonncau, Nucl.Phys. B 593, 398 (2001); M. Perez- Victoria, Phys. 
Rev. Lett. 83, 2518 (1999); M. Perez- Victoria, J. High. Energy Phys. 04, (2001) 032; O.A. Battistel and G. Dallabona, 
Nucl. Phys. B 610, 316 (2001); O.A. Battistel and G. Dallabona, J. Phys. G 28, L23 (2002); J. Phys. G 27, L53 (2001); A. 
P. B. Scarpelli, M. Sampaio, M. C. Nemes, and B. Hiller, Phys. Rev. D 64, 046013 (2001); T. Mariz, J.R. Nascimento, E. 
Passos, R.F. Ribeiro and F.A. Brito, J. High. Energy Phys. 0510 (2005) 019; J. R. Nascimento, E. Passos, A. Yu. Petrov, 
F. A. Brito, J. High. Energy Phys. 0706, (2007) 016; B. Altschul, Phys. Rev. D 70, 101701 (2004); A.P.B. Scarpelli, M. 
Sampaio, M.C. Nemes, B. Hiller, Eur. Phys. J. C 56, 571 (2008). 
[10] Q. G. Bailey and V. A. Kostelecky, Phys. Rev. D 70, 076006 (2004). 

[11] X. Xue and J. Wu, Eur. Phys. J. C 48, 257 (2006); H. Belich, M. M. Ferreira Jr, J.A. Helayel-Neto, M. T. D. Orlando, 
Phys. Rev. D 68, 025005 (2003); H. Belich, T. Costa-Soares, M.M. Ferreira Jr., J. A. Helayel-Neto, Eur. Phys. J. C 42, 
127 (2005); M. B. Cantcheff, Eur. Phys. J. C 46, 247 (2006); S.O. Vellozo, J.A. Helayel-Neto, A. W. Smith, and L.P.G. de 
Assis, Int. J. Theor. Phys. 47, 2934 (2008); Rodolfo Casana, Manoel M. Ferreira Jr and Carlos E. H. Santos, Phys. Rev. 
D 78, 025030 (2008). 

[12] Rodolfo Casana, Manoel M. Ferreira Jr and Carlos E. H. Santos, Phys. Rev. D 78, 105015 (2008). 

[13] Rodolfo Casana, Manoel M. Ferreira Jr, A. R. Gomes and Paulo R. D. Pinheiro, Eur. Phys. J. C 62, 573 (2009). 

[14] Rodolfo Casana, Manoel M. Ferreira Jr and J. S. Rodrigues, Phys. Rev. D 78, 125013 (2008). 

[15] J.M. Fonseca, A.H. Gomes, W.A.Moura-Melo, Phys. Lett. B 671, 280 (2009). 

[16] V. A. Kostelecky and CM. Pickering, Phys. Rev. Lett. 91, 031801 (2003); B. Altschul, Phys.Rev. D 70, 056005 (2004). 
[17] CD. Carone, M. Sher, and M. Vanderhaeghen, Phys. Rev. D 74, 077901 (2006); B. Altschul, Phys. Rev. D 79, 016004 
(2009). 

[18] M.A. Hohensee, R. Lehnert, D. F. Phillips, R. L. Walsworth, Phys. Rev. D 80, 036010(2009); M.A. Hohensee, R. Lehnert, 

D. F. Phillips, R. L. Walsworth, Phys. Rev. Lett. 102, 170402 (2009). 
[19] C. Adam and F. R. Klinkhamer, Nucl. Phys. B 657, 214 (2003). 

[20] F.R. Klinkhamer and M. Risse, Phys. Rev. D 77, 016002 (2008); F.R. Klinkhamer and M. Risse, Phys. Rev. D 77, 117901 
(A) (2008). 

[21] F.R. Klinkhamer and M. Schreck, Phys. Rev. D 78, 085026 (2008). 

[22] R. Montemayor and L.F. Urrutia, Phys. Rev. D 72, 045018 (2005); B. Altschul, Phys. Rev. D 72, 085003 (2005); B. 

Altschul, Phys. Rev. D 74, 083003 (2006); B. Altschul, Phys. Rev. D 72, 085003 (2005). 
[23] V. A. Kostelecky and M. Mewes, Phys. Rev. D 80, 015020 (2009). 

[24] A. P. Baeta Scarpelli and J. A. Helayel-Neto, Phys. Rev. D 73, 105020 (2006); N.M. Barraz, Jr., J.M. Fonseca, W.A. 
Moura-Melo, and J.A. Helayel-Neto, Phys. Rev. D76, 027701 (2007); H. Belich , J.L. Boldo, L.P Colatto, J.A. Helayel- 
Neto, A.L.M.A. Nogueira, Phys.Rev. D 68, 065030 (2003); G. de Berredo-Peixoto and I.L. Shapiro, Phys. Lett. B 642, 
153 (2006); A.P. Baeta Scarpelli, H. Belich, J.L. Boldo, L.P. Colatto, J.A. Helayel-Neto, A.L.M.A. Nogueira, Nucl. Phys. 
Proc. Suppl.127, 105 (2004); J.W. Moffat, Int. J. Mod. Phys. D 12 1279 (2003); F. W. Stecker and S.T. Scully, Astropart. 
Phys. 23, 203 (2005); E. O. Iltan, Eur. Phys. J. C 40, 269 (2005); Mod. Phys. Lett. A19, 327 (2004); J. High. Energy 
Phys. 0306 (2003) 016; O. Bertolami and D.F. Mota, Phys. Lett. B 455, 96 (1999); 

[25] M. B. Cantcheff, C.F.L. Godinho, A.P. Baeta Scarpelli, J.A. Helayel-Neto, Phys. Rev. D 68, 065025 (2003); H. Belich, T. 
Costa-Soares, J.A. Helayel-Neto M.T.D. Orlando, R.C Paschoal, Phys. Lett. A 370, 126 (2007); H. Belich, L.P. Colatto, 
T. Costa-Soares, J.A. Helayel-Neto, M.T.D. Orlando, Eur. Phys. J. C 62, 425 (2009); F.A. Brito, J.R. Nascimento, E. 
Passos, A.Yu. Petrov, Phys. Lett. B 664, 112 (2008); M. Gomes, T. Mariz, J.R. Nascimento, A.J. da Silva, Phys. Rev. D 
77, 105002 (2008); F.A. Brito, L.S. Grigorio, M.S. Guimaraes, E. Passos, C. Wotzasek, Phys. Rev. D 78, 125023 (2008). 

[26] T. Jacobson, S. Liberati, and D. Mattingly, Ann. Phys. 321, 150 (2006); T. Jacobson, S. Liberati, D. Mattingly and F.W. 
Stecker, Phys. Rev. Lett. 93, 021101 (2004); T. Jacobson, S. Liberati, and D. Mattingly, Nature 424, 1019 (2003). 

[27] A. Kobakhidze and B.H.J. McKellar, Phys. Rev. D 76, 093004 (2007). 

[28] J. I. Kapusta, Finite- Temperature Field Theory, Cambridge University Press, Cambridge, 1989 ; M. Le Bellac, Thermal 
Field Theory, Cambridge University Press, Cambridge, 1996. 



